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Abstract 

QCD sum rules for the nucleon are considered in complex q 2 plane 
with inclusion of the radiative corrections of the order a s . It is shown 
that the radiative corrections affect mainly the residue A 2 of the nu- 
cleon pole. Their influence on the value of the nucleon mass is much 
smaller. Following the ideas of Ioffe and Zyablyuk we expand the 
analysis to complex values of q 2 . This provides a more stable solu- 
tion. Varying the weights of the contributions of different dimensions 
by changing the value of the angle in the complex plane we find the 
value of the six-quark condensate which insures the best consistency 
of the right hand sides and left hand sides of the sum rules. The corre- 
sponding value of the six-quark condensate appears to be only about 
10% smaller then the one, provided by the factorization approxima- 
tion. The value of the four-quark condensate also appears to be close 
to the one, corresponding to the factorization assumption. The role 
of the gluon condensate and its possible values are discussed. 



I. INTRODUCTION 



The QCD sum rules (SR) method invented by Shifman et al. [I] succeeded 
in expressing the static properties of the hadrons in vacuum in terms of the 
expectation values of QCD operators. The SR approach was employed by 
Ioffe et al. dj in the case of nucleons. The basic point is the dispersion 
relation for the polarization operator 

11(g) = i I d A xe^ (0\Tj(x)j(0)\0) (1) 

of the local operator j(x) carrying the quantum numbers of the nucleon. The 
dispersion relation is considered at large values of \q 2 \, while q 2 < 0, where it 
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can be presented as a power series of q~ 2 (with the QCD condensates being 
the coefficients of the expansion). This presentation is known as the operator 
product expansion (OPE) At these values of q 2 the polarization operator 
can be presented also as a power series of the QCD coupling constant a s . 

Several lowest order OPE terms have been obtained in [21 E] ■ The leading 
contribution of the order q 4 In q 2 comes from the free three-quark loop. The 
higher terms contain the expectation values 

(0|w|0), (0\^G%G%\0), (0\qqqq\0), 

etc., with q and G a standing for the quark operators and for the gluon 
field tensor. The analysis of |3j included also the most important radiative 
corrections, in which the coupling constant a s is enhanced by the "large 
logarithm" lng 2 . The corrections (a s lnq 2 ) n have been included to all orders 
for the leading OPE terms [HUE]. This approach provided good results for 
the nucleon mass [3] and for the other characteristics of nucleons jfr. 

There are several traditional points of SR analysis. Following 3J one 
writes dispersion relations 

(i = q, I) for the components of the polarization operator 

n(?) = q^W(q 2 )+m\q 2 ) (3) 

with 7^ and / standing for the Dirac and unit matrices. The result of OPE 
on the left-hand side (LHS) of Eq. (2), is equalled to the hadronic contribu- 
tion on the right hand side (RHS), which is approximated by the standard 
"pole+continuum" model. The spectral density Imr?(A; 2 ) is approximated 
by the standard "pole+continuum" model 

- Im Ii\k 2 ) = X 2 N 5(k 2 - m 2 ) 

+ I^ 2 -^ 2 )Imn lOPE (A; 2 ). (4) 

7T 

Thus, the position of the lowest laying pole m, the residue X%, and the con- 
tinuum threshold W 2 are the unknowns of the SR equations. The standard 
Borel transform is usually carried out. 

However, inclusion of the lowest order radiative corrections beyond the 
logarithmic approximation made the situation somewhat more complicated. 
A numerically large coefficient of the lowest radiative correction to the leading 
OPE of the polarization operator (1) was obtained in A more consistent 
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calculation jH] provided this coefficient to be about 6. Thus, the radiative 
correction increases this term by about 50% at \q 2 \ ~ 1 GeV 2 , which are 
actual for the SR analysis. This "uncomfortably large" correction is often 
claimed as the most weak point of the SR approach 

The radiative corrections of the order a s In q 2 and a s for the contributions 
up to q~ 2 have been calculated by Ovchinnikov et al. [10] . These results were 
used for the calculation of the parameters m, X 2 and W 2 by finite energy sum 
rules method (FESR) [TTJ. Later the influence of the radiation corrections 
on the values of the nucleon mass was studied by Borel SR technique at fixed 
values of W 2 and A 2 , taken from FESR calculations [12] • Recently the QCD 
SR with the radiative corrections included have been used for determination 
of the value of quark scalar condensate However, the analysis of the 

role of the radiative corrections performed totally in framework of the Borel 
transformed SR have not been done until now. This is done in the present 
paper. We investigate also the sensitivity of the nucleon characteristics to 
the variation of the values of the condensates of higher dimensions. 

These two points are connected with each other. While the condensate 
(0|gg|0) of dimension d = 3 is known with good accuracy, there is no way to 
obtain the four- and six-quark condensates (d = 6,9) from an experiment. 
The calculation of these expectation values requires some additional assump- 
tions. As it stands now, the only approach used in the calculations is the 
factorization hypothesis, which assumed the domination of the vacuum states 
in the sum over the intermediate states pQ. In this approximation expecta- 
tion values of the products of four and six quark operators are expressed in 
terms of the scalar operator (0|gg|0). However, the correction of the order 
a s ln(q 2 / fi 2 ) to this condensate differs from the sum of such corrections to 
the condensates (0|gg|0) [TUj. Thus, the factorization hypothesis should be 
combined with the definition of the normalization point (scale) at which it 
is assumed to be true. 

We study also the dependence of the nucleon parameters on the value of 
the gluon condensate (0|— G^G^O) (d = 4). There are indications that 
this matrix element can be smaller [2] or larger than the standard value 
obtained in [To] . 

Analysis of the SR at complex values of q 2 = \q 2 \e tip was employed first 
by Ioffe and Zyablyuk for investigation of hadronic r-decay. Varying 
the value of angle ip in complex plane of q 2 one can change the weights of 
different OPE contributions. Some of the terms can be eliminated in such 
a way. By using this technique one can expect to obtain more reliable and 
stable results. The hadron parameters should not depend on ip. We apply 
this approach for the nucleon channel. 

We recall the original form of QCD SR in Sect. II. Corrections of the 
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order a s are included in Sect. III. The corresponding equations in complex 
q 2 plane and the analysis are presented in Sect. IV. 



II. QCD SUM RULES AT REAL VALUES OF q 2 



For the specific calculations we use the "current" j, which enters Eq. (1) 
in the form [2] 



with u and d standing for the quark fields. C is the charge conjugation 
operator, while a, b, c are the color indices. The lowest OPE terms of the 
operators U q (q 2 ) and IT 7 (q 2 ) introduced by Eq. (3) can be presented as 



Here the lower indices show the dimensions of the condensates, contained in 
the corresponding terms, Aq is the contribution of the free quark loop. 

The leading OPE terms Aq, A^ and B% contain divergent integrals. It 
is sufficient [2] to carry out regularization in the simplest way, i.e., just by 
introducing an ultraviolet cutoff C u . Direct calculation with all the radiative 
corrections being neglected provides j2] 



with the standard notation G 2 = G a G a , while Q 2 = -q 2 > 0. 

The condensates of the lowest dimensions d = 3,4, i.e., (0|gg|0) and 
(0|— G 2 |0) can be obtained from the Gell-Mann-Oakes-Renner relation [18J 
and from the meson sum rules ^E]. The higher OPE terms Aq and B$ 
contain the expectation values of the products of four and six quark operators. 
There are neither experimental no rigorous theoretical data on the four- 
quark condensates of the general form (0\qTiqqT 2 q\0), with r 12 acting on 
Lorentz and color indices. The same refers to the six-quark condensates. 
The standard approach is the factorization approximatiom pp expressed by 
equation 



jW = 



u a (x)C^u b (x) 75 ^d c ( 



(5) 



A + A 4 + A 6 + A 8 , II 7 = B 3 + B 7 + B 9 . 



(6) 




<0|gr ig gr 2 g|0> = 

= tf- 2 [(T!rr 1 '&r 2 ) -Tr^r^KOlgglO) 2 



4 



with q standing for u and d quarks and N = 12. This equation presents all 
the four-quark condensates in terms of the scalar expectation values (0|gg|0). 
However all Lorentz structure with Ti = I^ contribute due to the second term 
in brackets on the RHS of the latter equation (several examples are presented 
in £Q). Similar relation can be written for the six-quark condensate. Also, the 
factorization approximation for the quark-gluon condensate, which enteres 
the expression for B 7 provides 

(0|g^G 2 g|0> = (0|g#> (0| - G 2 |0> . 

7T 71 



Thus, the higher order terms are 

2 
3Q 



^6 = 777^2 



A 8 = -^/i 2 ((0|gg|0>) 2 , 

5r = T ^(0|gg|0)(0|^G 2 |0), 

B 9 = 32 7 r 2 ^i I ((0|gg|0>) 3 . (8) 

Here the terms A e , B 7 and B 9 are presented in the factorization approxi- 
mation. In A s we have taken into account that the sum rules for the baryon 
resonances provide JTHj 

(OlM-G^^lO) = /, 2 ((0|gg|0» 2 (9) 

71 r Z 

with /i 2 « 0.8 GeV 2 . 

In the analysis, carried out in the most important radiative corrections 
of the order a s In q 2 have been included. These contributions were summed 
to all orders of (a s lng 2 ) n , being expressed in terms of the factor 



' lnQ 2 /A 2 \ 
^ In /i 2 /A 2 J 



4/9 



(10) 



with A w 150 MeV is the QCD scale, while \i is the normalization point. 
It was assumed in jH] that the physical characteristics are normalized at 
fi = 500 MeV. The power 4/9 reflects the so-called anomalous dimension. 

The next traditional point is the Borel transform 

Q = — q — > 00, n — > 00, — = M , 

n 
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converting a function of q 2 into the function of M 2 . It removes the terms 
depending on C u in the expressions provided by Eq. (7), and also makes 
the "pole+continuum" model more reliable, suppressing the contributions 
of the heavier states by the factor e~ k l M . Actually, in order to deal with 
the values of the order of unity (in GeV units), it is convenient to use the 
operator B* = 32n 4 B. 

After these manipulations the Borel transformed SR are [3] 

C\M 2 ) = R q (M 2 ), L X {M 2 ) = R\M 2 ) (12) 

with 

C q = A + I 4 + Ae + A 8 , L 1 = B 3 + B 7 + B 9 . (13) 

Here, as well as in Eq. (6), the lower indices show the dimensions of the 
condensates. The terms on the RHS of Eqs. (13) are 

~ M G E 2 ~ bM 2 E ~ 4 
A) = L , At = 4L , A e = - a L, 

i 8 = -^« 2 , B 3 = 2aM 4 E 1 , 

Br = --, B9 = — -4 (14) 
12 81 7r M 2 K J 

with the notations E { = Ei(W 2 /M 2 ),i = 0, 1, 2 and 



O.s 



a = -(2vr) 2 (0|gg|0) and b= (2tt)^<0|— G 2 |0). (15) 
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The conventional values of these condensates at fi = 0.5 GeV are a = 0.55 GeV 3 
and b = 0.50 GeV 4 . The functions 



E (x) = 1 - e~ x , E^x) = 1 - (1 + x)e 



■X 



E 2 (x) = 1- (^l + x + ^Pje-* (16) 

account for the contribution of the continuum, Ei(x) — - > 1 for W 2 — ► oo. 
Thus, the RHS of Eqs. (12) 

R\M 2 ) = X 2 e- m2 / M \ R\M 2 ) = mA 2 e" m2/M2 (17) 

describe only the contributions of the nucleon pole with A 2 = 327r 4 A 2 v . 

The LHS of Eq. (12), being calculated as OPE series, becomes the bet- 
ter approximation, while the value of M 2 increases. On the other hand, 
"pole+continuum" model for the RHS becomes increasingly accurate for 
smaller values of M 2 . The interval 

0.8 GeV 2 < M 2 < 1.4 GeV 2 , (18) 



6 



where both approximations are expected to be true was found in The 
values of the parameters m, A 2 and W 2 have been obtained by minimization 
of the functional 



h{M 2 ) 



i=q,I \ 



C\M 2 ) - R\M 2 ) 
£*(M 2 ) 



(19) 



in the interval (18), providing 



m = 0.931 GeV , A 2 = 1.86 GeV 6 , 
W 2 = 2.09 GeV 2 . 



(20) 



Account of the anomalous dimensions is not very important. Putting L = 1 
one can find 



Note that straightforward using of this approach can provide mislead- 
ing results. The solution of (19) corresponding to minimization appears to 
be unstable with respect to the values of the QCD condensates. Moreover, 
sometimes the functional (19) may have few different minima. This is illus- 
trated by Fig. 1, where the solution corresponding to the second minimum is 
shown by the dashed line. For the values of b which do not exceed strongly 
the conventional value b = 0.50 GeV 4 the best \ 2 fitting is provided by the 
first solution shown by the solid line. However, for larger values of b the min- 
imization procedure makes us to jump to the second solution with smaller 
value of the nucleon mass m = 0.6 GeV. Consider, for example, the SR with 
somewhat larger value of the gluon condensate. At b = 0.65 (which is 30% 
larger than the conventional value) minimization of fi (19) provides much 
smaller values m ~ 0.6 GeV and W 2 = 1 GeV 2 -see Fig. 1. Such solution can- 
not be treated as a physical one because of the small value of W 2 . Indeed, 
the contribution of the continuum (treated approximately) exceeds the con- 
tribution of the pole, which is treated exactly. Thus the "pole+continuum" 
model for the RHS of the sum rules has no sense. Although the accuracy 
of the physical solution does not change much with b, it will not be noticed 
by x 2 minimization procedure since unphysical solution is even more accu- 
rate — see Fig. 1. The unphysical solution with small values of m and W 2 
can be traced by successive inclusion of the condensates of higher dimension. 
Including only the condensates of dimensions d = 3,4 we find a trivial solu- 
tion m = 0, A 2 = 0, W 2 = 0. Inclusion of the condensate with d = 6 still 
keeps m = 0,W 2 = 0, but A 2 = |a 2 = 0.4 GeV 6 . Inclusion of the higher 
condensates provides small nonzero values of m and W 2 . 



m = 0.930 GeV, A 2 = 1.79 GeV 6 , 
W 2 = 2.00 GeV 2 . 



(21) 
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To avoid this situation we add analysis of the functions 



m\(M 2 ) 
m 2 2 {M 2 ) 



M 4 dC q (M 2 



£i(M 2 ) dM 2 

M 4 dC^M 2 ) 
£ 7 (M 2 ) dM 2 



ms (M=) = (22) 

Note that these functions depend also on the threshold value W 2 . For the 
solution of the SR equations (12) it should be 

m 2 (M 2 ) = m 2 2 (M 2 ) = m 2 (M 2 ) = const = m 2 . (23) 

Note different meaning of the masses my and m 3 . While my determine the 
position of the lowest pole, which approximates the q 2 dependence of (3), the 
difference m 3 — m reflects rather the admixture of the negative parity state 
l/2~ to the nucleon. 

III. CORRECTIONS OF THE ORDER a s AT REAL VALUES OF q 2 

Inclusion of the corrections of the order a s to the polarization operator 
modifies the expressions on the RHS of Eqs. (7), (8) to 



64tt 4 /i 2 



= --— g 4 in^- 



/ , Ua s la, Q 2 \ 

x 1 H In — , 

\ 12 tt 2 vr /J 2 ) ' 

2(0|™|0)/ 5« s la, Q 2 
Ae " 3 Q 2 7 

4tt 2 /i 2 V 2 7T 



S 3 = -^T^Q 2 ln^(l + ^), (24) 



obtained in 0110] (see Appendix A). The radiative corrections to the other 
terms are not included since the values of the corresponding condensates are 
known with poor accuracy. Strictly speaking this is true for the four-quark 
condensate as well, which is obtained in the factorization approximation only. 
However, we shall assume the latter hypothesis for the four-quark condensate, 
advocated recently in |20| . Note that the correction of the type a s ln(g 2 //x 2 ) 
to the term B-i vanishes due to the cancellation of the correction ^— In % 
to the condensate (0|gg|0) and that coming from the free quark system. If 
the factorization approximation would have been true for all values of /i 2 , we 
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would find +— In % for the last term in brackets of the formula for Ar. One 
can see that the value of the coefficient is another one. Thus, the factorization 
assumption cannot be true for all values of fi 2 . 1 

In Ref. jinj these formulas have been used for investigation of the nucleon 
parameters with the finite energy sum rules technique [TT] . It was shown that 
inclusion of radiative correction leads to a noticeable reduction of the nucleon 
mass. Similar behavior was obtained in the analysis carried out in ^2] i n 
framework of the Borel transformed sum rules at fixed values of W 2 . 

Using (24), we find for the corresponding contributions to the Borel trans- 
formed SR (see Appendix B) 



A (M 2 ,W 2 ) = M G E 2 



a* 



1 + — — -hi — 



7T 
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M 4 W 2 



' 3W 2 \ 
+ AM 2 ) 




+ M 6 £[-W 2 /M 



A 6 (M 2 ,W 



B 3 (M 2 ,W 



with Ei = Ei(W 2 /M 2 ) 




n=l 



n ■ m 



(25) 



Expression for B3 coincides with that presented in ^2j- Two other ones 
are more general than those of ^2], being true for any value of fi 2 , and 
accounting for the contribution of continuum in the main radiative correction 
to the term of the sixth dimension. Besides that some terms in the expression 
for Aq differ from those in ^2] • The numerical difference is, however, not very 
important. 

One can see that the free loop term Aq suffers the largest radiation cor- 
rections. If the threshold W 2 is large enough, Aq provides the leading contri- 

1 Since the value of ln(M 2 //Lt 2 ) is not large and we are studying the role of the whole 
0{a s ) corrections, and not only the part (given by the renormgroup) which is responsible 
for the scale dependence of the polarization operator, here we keep the logarithmic a s In Q 2 
contribution but, to avoid the double counting, neglect the anomalous dimension factor 
Li. 
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bution to the structure C q . Hence, neglecting the influence of radiative cor- 
rections on the value of W 2 , one can see that the denominator of the function 
m^(M 2 ) (22) increases due to the radiative corrections. Thus, in agreement 
with ^OJ^], the latter diminish the value of nucleon mass. However, the 
present analysis which includes possible modification of all the parameters 
(m, A 2 , W 2 ) by the radiative corrections, provides somewhat different result. 

To carry out quantitative analysis, we need to clarify the argument of the 
running coupling constant 

An 

°-<*> = Impm 5 ) • (26) 

In PHI the radiative corrections have been obtained at a s = const. Since the 
momenta in the loops corresponding to the radiative corrections are of the 
order of q, it is reasonable to assume a s = a s (Q 2 ) in Eq. (24). After the Borel 
transform we obtain a s = a s (M 2 ) in Eq. (25). Since we are considering M 2 
of the order 1 GeV , we can put a s = a s (1 GeV 2 ) = 0.35. Minimization of 
the functional (19) provides 

m = 0.94 GeV , A 2 = 2.00 GeV 6 , 

W = 1.90 GeV 2 , (27) 

while for a s = a s (M 2 ) 

m = 0.94 GeV, A 2 = 2.11 GeV 6 , 

W 2 = 2.00 GeV 2 . (28) 

In Fig. 2 we present the values of the parameters as the functions of 
a s . One can see that the radiative corrections affect mostly the value of 
the residue A 2 . The consistency of the LHS and RHS of the SR is shown in 
Fig. 3(a). 

Note that in [Til] the SR have been considered at fixed value W 2 = 
2.5 GeV 2 , providing somewhat larger residue value A 2 ~ 3 GeV 2 . Assuming 
W 2 = 2.5 GeV 2 in Eqs. (12), we would move to larger values A 2 = 2.6GeV 2 
in our solution. This is close to the result of , although our procedures of 
inclusion of radiative corrections differ in several points. 

We can look for the solution in which the discrepancy between the val- 
ues of rrii(M 2 ) defined by Eq. (22) is minimized. It can be obtained by 
minimization of the functional 

f 2 {M 2 ) = h{M 2 ) + f rn {M 2 ) (29) 

with 

n ^ ( m iM 2 ) -m) 2 

UM 2 ) = E 1 1 ^ L , (30) 
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while f\(M 2 ) is determined by (19). For a s = the procedure provides 
m = 0.93 GeV, A 2 = 1.82 GeV 6 , W 2 = 2.02 GeV 2 , while for a s = 0.35 we 
find 

m = 0.94 GeV , A 2 = 1.98 GeV 6 , 

W 2 = 1.88 GeV 2 . (31) 

Thus, minimization of the functionals fi(M 2 ) (19) and f 2 (M 2 ) (29) leads 
to close results. The functions mj(M 2 ) are presented in Fig. 3(b). 

IV. SUM RULES IN THE COMPLEX PLANE OF Q 2 

Now we consider the SR at complex values of 

Q 2 = S 2 e i(f (32) 

with the real values of S 2 > and — n < ip < n. The Borel transform (11) 
will be carried out with respect to S 2 . It is reasonable to consider the SR for 
the real part of the operators IP(g 2 ) (3). 

A. Possible values of ip 

Note first that our "pole+continuum" model makes sense only for cos tp > 
0, i.e., — | < </? < |. This becomes clear while considering the Borel trans- 
form of contribution of the nucleon pole to the RHS of Eq. (2) 

B * f ^y 2 - m 2 ) ^ = A a * (33) 

b J m 2 + S 2 e^ v ' 



,2 



with 9 — ^2 sin (p. Thus the contributions of heavier states are suppressed 
with respect to the lowest one only if cos(/? > 0. 

Also, the values of <p should not be too close to ±7r/2. This is because the 
leading OPE term B 3 obtains the factor cosy? in the complex plane, while 
the higher term B 9 is multiplied by cos2y?. Hence, the convergence of OPE 
series becomes worse at (p close to ±ir/2. That's why we focus on the values 

< V < \ ■ (34) 

Analysis at negative values of <p will not provide new data, since all the 
functions involved are the even functions of ip. 

B. Basic equations in the complex plane 
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Following the previous discussion, we must present the Borel transformed 
dispersion relations for the real parts of the operators U q (q 2 ) and Il / (g 2 ). For 
the Borel transforms of the contributions (7), (8) we find 

M 6 e 2i,p - M 2 a 

BA « = ^ ^ = ^<0|^|0>, 

BA 6 = ^e-^(0\qqqq\0), 

BA 8 = g=±j e-**(0\m* G% y <v<z|0) , 

BB s = ~e^(Q\gq\0) , 
BB 7 = ^(0\q^G 2 q\0) , 

24 7T 

272e~ 2iv a 

BB 9 = glM2 (0\qqqqqq\0) ■ (35) 

The contributions A , A A and B 3 contain the quark loops, thus requiring 
inclusion of the contribution of the continuum. The latter can be evaluated 
as continuum, which takes the form 

_ 7 dk 2 (k 2 r _ f dk 2 c {k 2 c ) n e^ 

n J k 2 + S 2 e^ ~ J k 2 + S 2 [ ' 

w 2 wt 

with n — 0, 1, 2, and 

W 2 = W 2 e~^. (37) 
The corresponding Borel transforms are thus 

BR n = e m *M 2 (l - E n {jpj^ . (38) 

Using these results, and including the radiative corrections (25) we can 
present the SR in the form 

Re = Rl, Re = R J C (39) 

with 

Rl = \ 2 e~^ cos (^sin^, 



2 



Ri = m\ 2 e-^ cosip cos -^siny?), (40) 
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while 



C q c (M 2 , W 2 )=e 2 ^A (M 2 , W 2 ) + A A (M 2 , W 2 ) 

+ e- llp A 6 (M 2 ) + e- 2llp A 8 (M 2 ), (41) 

Ci(M 2 , W 2 ) = e iip B 3 (M 2 , W 2 ) 

+ e~ iip B 7 + e- 2ilp B 9 (M 2 ) (42) 

with W 2 defined by Eq. (37), while the contributions Ai and Bi are given by 
Eq. (25). The functions 

2 /,,2x M 4 e iip dC q (M 2 ) 

ml(M 2 ) = Re 



m 2 (M 2 ) = Re 



£i(M 2 ) dM 2 
M 4 e^ dCHM 2 ) 



2V ; C J {M 2 ) dM 2 ' 

Re£ 7 (M 2 ) 

m: ' iM » = rtt^ (43) 

generalize Eq. (22) for the complex plane, and Eq. (23) should have been 
true for the exact solution 



C. Dependence of hadron parameters on ip 

Now we try to find the nucleon parameters by minimization of the func- 
tional 

which is generalization of Eq. (19) for the complex plane of Q 2 . Since the 
M 2 interval of stability may vary with tp, we present the results for the two 
intervals, e.g., the one given by Eq. (18), and also 

1.0 GeV 2 < M 2 < 1.8 GeV 2 . (45) 

The results are shown in Fig. 4. One can see that the largest stability is 
obtained at ip < |. Also, the nucleon residue appears to be more sensitive 
to the value of tp, then the values of the nucleon mass and of the continuum 
threshold W 2 . 



D. Dependence on the value of six-quark condensate 

As we said earlier, the less reliable value of the condensates involved is 
that of the six-quark condensate, which determines the contribution B$ on 
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the RHS of Eq. (42). While all the results obtained above are obtained in 
factorization approximation, now we put 



(0|(gg) 3 |0) =((0|gg|0)) 3 r ?6(? , 



(46) 



with 7]Q q = 1 under the factorization assumption and check dependence of 
the parameters on ip at various values of r) 6g . The results are presented in 
Fig. 5. One can see that the value of the nucleon mass m exhibits larger 
stability with respect to changes of p> and rje q then the values of A 2 and W 2 . 



Note that at ip = n/4 the contribution Aq on the RHS of Eq. (41), corre- 
sponding to the free-quark loop also turns to zero. This is not necessary 
a weak point, since the system of three free quarks contributes rather to 
continuum then to the lowest laying pole. 

At if = 7r/4 the consistency of LHS and RHS of Eqs. (39) is not very good. 
It becomes much better at smaller values of y? in the domain of stability in 
(M 2 , ip) plane. 



Now we shall look for the solution which provides the best consistency of 
the LHS and RHS of Eqs. (39) in certain domains of the values of M 2 and ip. 
In Table I we present result obtained by minimization of the functional (44). 
One can see that the consistency ("% 2 per point") becomes much better if 
we limit the interval of the values of the angles to 



The results obtained in the interval (45) are close to those found for the 
interval (18). In this case "x 2 P er point" is somewhat smaller than in the 
latter one. However we consider the results obtained in traditional interval to 
be more reliable, since on the upper end of (45) the contribution of continuum 
is not suppressed stronger than that of the pole. Thus, the limits of stability 
of SR in (M 2 ,(p) plane are given by Eqs. (18) and (48) with the values of 
the parameters being 




E. Stability of solution in (M 2 ,ip) plane 



< (p < - . 



(48) 



m = 0.94 GeV , A 2 = 2.0 GeV 6 , 
W 2 = 1.9 GeV 2 . 



(49) 
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Consistency of LHS and RHS of Eqs. (39) is illustrated by Fig. 6 for the 
hadron parameters presented by Eq. (49) for several values of tp. In Fig. 7 
we show the functions rrii(M 2 ) defined by Eq. (43) for several values of tp 
and dependence of the values mj(M 2 ) on tp for several values of M 2 . This 
illustrates the weak dependence of the values of the hadron parameters on 
the values of M 2 and tp in this region. 

Note that similar results can be obtained by minimization of the func- 
tional 

f c 2 {M 2 ) = r x {M 2 ) + f m {M 2 ) (50) 

with f m and /f defined by Eqs. (30) and (44). These data are presented in 
Table II . The functions rrn(M 2 ) are shown in Fig. 8. 

F. Values of the condensates of high dimension 

Now we can try to use the SR for diminishing of the uncertainties of the 
values of the QCD condensates. We put 

(0|(gg) 2 |0) = ((0|gg|0)) 2 77 4g , 

(0\-G 2 \0) = (2tt)- 2 Mg (51) 

71 

with b = 0.50 GeV 4 , and try to determine the values of rj4 q , t)q, and r) 6q (46) 
from the SR. Deviations of the parameters 77 4(? and r] 6q from unity characterize 
the violation of factorization hypothesis. Possible deviation of parameter rja 
from unity is due to uncertainties in the value of the gluon condensate. The 
terms A±, A e and B 9 on the RHS of Eqs. (41) and (42) obtain the factors 
Vg, V^q an d V6q correspondingly. 

Now we try to solve SR equations treating the hadron parameters m, A 2 , W 2 
and also the values r)4 q , r]e q and r/c as the unknowns. 

Four-quark condensate. Here we put rj 6q — rja — 1 and try to determine 
the value of r/4 q . Minimization of (44) provides 

m = 0.964 GeV , A 2 = 2.11 GeV 6 , 

W 2 = 1.99 GeV 2 , r] 4q = 0.93. (52) 

Thus, the QCD sum rules prefer only small deviation from the results of the 
factorization approximation. 

Six-quark condensate. Now we try to determine the value of i]Q q , putting 
f}/± q = t]q = 1. Minimization of (44) provides 

m = 0.937 GeV , A 2 = 1.97 GeV 6 , 

W 2 = 1.90 GeV 2 , r] 6q = 0.90, (53) 
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demonstrating small deviations of the value r] eq from the factorization hy- 
pothesis value r) 6q — 1. 

Gluon condensate. Now we fix r] iq = r]Q q = 1 and look for the dependence 
of the solution on the value of the gluon condensate. The value of the nucleon 
mass appears to be rather sharp function of r/c- For example, at r/c = 2 we 
find m = 0.74 GeV. 

Note that the functions m^(M 2 ) (22) is much more sensitive to the value 
of the gluon condensate than the "pole masses" mi )2 (M 2 ). Variation of the 
parameter r/c in the interval < r/c < 2 changes the values of by about 
60MeV, while m 3 changes by a factor 1/2. 

Treating r\ G as an unknown parameter of the SR equations we find that 
minimization of (44) takes place at t]q = 2.2 and at the underestimated value 
of the nucleon mass m ~ 0.67 GeV. On the other hand, one needs rjo ~ 1 to 
obtain the value of the nucleon mass close to the observable one. 

V. SUMMARY 

We carried out the analysis of QCD sum rules for nucleons including 
the lowest order radiative corrections. This is the first analysis performed 
totally in framework of SR approach. We show that the radiative corrections 
modify mainly the values of the nucleon residue, while that of the nucleon 
mass suffers minor changes. 

Our analysis was carried out for the real values of q 2 and in the complex 
q 2 plane. We found the region of stability in (M 2 , ip) plane in the latter case. 
Our main result is expressed by Eq. (49), being illustrated by Fig. 6. 

We used the nucleon sum rules to clarify the values of the condensates 
of high dimension. We showed that the sum rules require the deviations 
of the four-quark condensates from the factorization value to be very small. 
Similar deviations of the six-quark condensates do not exceed 10%. As to the 
uncertainties of the gluon condensate, the values, which exceed the standard 
one by 30% provide the value of the nucleon mass m < 0.9 GeV. The greater 
values of the gluon condensate may cause problems in description of the 
nucleon. A more detailed analysis of the limitations on the gluon condensate 
value, coming from the nucleon sum rules will be published elsewhere. 
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APPENDIX A 



Here we describe the renormalization procedure, which leads to Eqs. (24). 
The point have not been considered in detail in previous publications. Until 
we do not include the radiative corrections to the quark loops (terms A , 
A 4 and B 3 in Eqs. (7), (8)), renormalization is not important, since the 
contribution Q 4 ln(C 2 ) is eliminated by the Borel transform. In other words, 
the Borel transform carries out the renormalization automatically. However, 
the situation becomes more complicated if the corrections of the order a s lnQ 2 
are included. 

Consider first the three-quark loop. While the radiative corrections are 
neglected, we can carry out the renormalization of the function Aq(Q 2 ) by 
subtracting the two lowest terms of the Tailor expansion at a normalization 
point Q 2 = fi 2 . This provides for the renormalized contribution 

(we omitted the terms, which will be killed by the Borel transform). The 
contributions which include the radiative corrections to Eq. (Al) are renor- 
malized in the same way. This leads to the equality, presented by the upper 
line of Eq. (24). 

Similar procedure can be applied to the calculation of the renormalized 
contributions Aq and S3. The result is presented by Eq. (24). There are no 
loops in the case of Aq. Thus, only the radiative corrections are renormalized. 
There are no corrections of the order a s lnQ 2 in the case of B 3 . Hence, 
only the two quark loop was renormalized. Unlike the case of A , the 
renormalization was not important here due to the Borel transform. 

APPENDIX B 

Here we calculate the contribution of the term A determined by Eq. (24) 
to the LHS of Eq. (12). This means that we must calculate the Borel trans- 
form, subtracting the contribution of the continuum. 



The renormilized contribution is 

'64vr 4 ~* *~ [i 2 



x [l + Ii^-I^ln^]. (Bl) 

1 12 71 2 71 fi 2 I 



Employing [2J we find 

S(Q 4 lnQ 2 ) = -2M 6 
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B(Q A In 2 Q 2 ) = -AM 2 (In M 2 - lE + - 
with 7s 0.577 being the Euler constant. Thus 

71 a. 



(52) 



5* A 



M 6 I 1 



Or 



„ 6 a s / M 2 
- M — In 

7T \ 



12 7T 



(S3) 



with the last term of the RHS originating from the last term in brackets on 
the RHS of Eq. (Bl). Recall that B* = 32n 4 B. 

Following Eq. (4) we find for the contribution to the RHS of Eq. (13) for 



MM 2 ,™*) = M«E 2 [1+ 7 ^ 

12 IT 



lE + -\-T(M 2 ,W 2 )\, 



with 



T{M\W 2 ) = - [ dk 2 k A \n^e~ k2 ' M2 (55) 
w 2 

describing the contribution of continuum, corresponding to the last term on 
the RHS of Eq. (B3). 

Introducing 



Z ~ M 2 ' 
and denoting k 2 = W 2 x, we evaluate 

T(M 2 ,W 2 ) = M 6 



(56) 



x 



z 2 \ W 2 
1 + z + — J e~ z ln — + t(z) 



with 



t(z) 



/ 



dxx \nxe~ 



{Bl) 



(58) 



Employing integration by parts three times, we obtain 

3 + z _ y 



«w =/ 



dx- 



x 



+ 



= -Ei(-z) + ——e 



(B9) 
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with 

Ei(x) =V.p. [ dy— = Y ^— + 7^ + 111 x (BIO) 

being the integral exponential function [21]. Combining Eqs. (B4), (B7) and 
(B9), we come to Eq. (25) of the text. 
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Figure captions 

FIG. 1. Dependence of solutions of Eq. (12) on the value of the gluon con- 
densate. The solid lines correspond to 6=0.50, the dashed lines are for the 
case 6=0.65. The nucleon parameters are related to those, determined by 
Eq. (21), i.e., m = 0.93 GeV, A 2 , = 1.79 GeV 6 , W 2 = 2.00 GeV 2 . 

FIG. 2. Dependence of ratios m/m , A 2 /Aq, and W 2 /Wq on the values of a s . 
Long-dashed line shows "x 2 per point". Parameters m , A 2 , and Wq are the 
same as in Fig. 1. 

FIG. 3. (a) Consistency of LHS and RHS of Eq. (12). Hadron parameters are 
given by Eq. (28). Solid and dashed lines show the ratios ^[^2] and ^rpfsy, 
correspondingly, (b) The functions rrii(M 2 ), i=l,2,3, defined by Eq. (22). 

FIG. 4. Dependence of the hadron parameters on the value of (p. Figures (a) 
and (b) correspond to the values of Borel masses in the intervals (18) and 
(45). Parameters m , A 2 , and Wq are the same as in Fig. I. 

FIG. 5. Dependence of the hadron parameters m (a), A 2 (b) and W 2 (c) on 
the value of ip for several values of r\§ q . The latter are shown by the numbers 
on the lines. 

FIG. 6. Consistency of LHS and RHS of Eq. (39) with the hadronic param- 
eters presented by Eq. (49) for <p = (solid line), ip = (dashed line) and 
(p = | (dotted line). 

FIG. 7. (a) The functions mj(M 2 ), i— 1,2,3, for ip—0 (solid line) and ip = | 
(dashed line). The functional f\(M 2 ) (19) is minimized, (b) Dependence of 
the functions m^M 2 ) on the value of <p for M 2 =0.8 GeV 2 (solid line) and 
for M 2 =1.4 GeV 2 (dashed line). 

FIG. 8. The functions m^M 2 ), i=l,2,3, for <p=0 (solid line) and ip = f 
(dashed line). The functional /|(M 2 ) (50) is minimized. 
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TABLE I. Nucleoli parameters obtained by minimization 
of the function f£ defined by Eq. (44). 



M 2 (GeV 2 ) 





x 2 


m (GeV) 


A 2 (GeV 6 ) 


W 2 (GeV 2 ) 


0.8 - 1.4 


0-f 


6.5 


0.93 


1.88 


1.84 


1.0 - 1.8 


o-f 


2.0 


0.96 


1.98 


1.86 


0.8 - 1.4 


°-f 


0.8 


0.94 


2.02 


1.93 


1.0 - 1.8 


o-f 


0.4 


0.96 


2.00 


1.88 



TABLE II. Nucleon parameters obtained by minimization 
of the function /| defined by Eq. (50). 



M 2 (GeV 2 ) 


<J 




x 2 


m (GeV) 


A 2 (GeV 6 ) 


W 2 (GeV 2 ) 


0.8 - 1.4 


0- 


7T 

4 


7.9 


0.94 


1.93 


1.88 


1.0 - 1.8 


0- 


7T 

4 


3.7 


0.97 


1.99 


1.87 


0.8 - 1.4 


0- 


7T 

8 


2.5 


0.94 


2.04 


1.95 


1.0 - 1.8 


0- 


7T 

8 


2.7 


0.96 


2.02 


1.90 
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